Universal shapes formed by two interacting cracks 



O 



in 

(N 



o 

s 
a 

O 

o 



> 

<N 

o 
o 



^ 



Melissa L. Fender,^ Frederic Lechenault,^'^ and Karen E. Daniels^ 

^Department of Physics, NC State University, Raleigh, NC, USA, 27695 

^LCVN, UMR 5587 CNRS-UM2, Universite Montpellier II, place Eugene Bataillon, 34095 Montpellier, France 

(Dated: 25 May 2010) 

We investigate the origins of the widely-observed "en passant" crack pattern which forms through interactions 
between two approaching cracks. A rectangular elastic plate is notched on each long side and then subjected 
to quasistatic uniaxial strain from the short side. The two cracks propagate along approximately straight paths 
until they pass each other, after which they curve and release a lenticular fragment. We find that for materials 
with diverse mechanical properties, the shape of this fragment has an aspect ratio of 2:1, with the length scale 
set by the initial crack offset s and the time scale set by the ratio of s to the pulling velocity. The cracks have a 
universal square root shape which we understand using a simple geometric model of the crack-crack interaction. 

PACS numbers: 62.20.mm, 46.50.+a, 



Brittle failure through multiple cracks occurs in a wide va- 
riety of contexts, from microscopic failures in dental enamel 
|1 1 and cleaved silicon \li to geological faults ISj |4j| and 
planetary ice crusts O. In each of these situations, with 
complicated curvature and stress geometries, pairwise in- 
teractions between approaching cracks nonetheless produce 
characteristically-curved fracture paths known in the geologic 
literature as en passant cracks (61. While the fragmentation of 
solids via many interacting cracks has seen wide investigation 
[7 - 15 1, less attention has been paid to the details of individ- 
ual crack-crack interactions. Despite extensive observations 
of the phenomenon of overlapping cracks in geologic settings 
El mo [161 as well as controlled laboratory experiments ifTvl - 
[22l and simulations 123^251 , an understanding of how similar 
shapes arise from different dynamics remains lacking. 

Through experiments and a simple geometric model, we 
elucidate the mechanism behind these two-crack interactions 
and quantify the scaling properties which allow for the same 
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FIG. 1 : (a) Experiment schematic, with Side 1 stationary and Side 2 
pulled at constant velocity v. The point at which the two crack paths 
first cross their mutual perpendicular sets t = 0, and [x^y] = 0. (b- 
f) Detail of gelatin sheet showing dynamics. Detail of photoelastic 
response of gelatin sheet (g) before and (h) during crack curvature. 



shape in diverse situations. We choose a geometry for which 
two-crack interactions are the dominant failure mode of the 
material: a uniaxial strain test on two prepared notches on 
opposite lateral sides of an elastic sample. As the ends of 
the sample are pulled apart quasistatically, the two notches 
propagate inward as single cracks until they pass and begin 
to curve towards each other, cutting out a lenticular fragment 
(see Fig.[T]3-f). The direction of this curvature is controlled by 
the ratio of mode I and mode II stress intensity factors Iil6l . 

Using materials varying in bulk modulus, Poisson ratio, and 
heterogeneity, we uncover a universal, scale-invariant shape 
and dynamics for this two-crack interaction. The shape of the 
final fragment exhibits rate-, material- and scale-invariance, 
with an aspect ratio (length: width) of F = 2 and an approxi- 
mately square root shape which begins at the point where the 
two cracks pass each other during propagation. We observe 
that the length scale of the fracture curvature is set by only 
the initial notch offset s, and that the time scale is set by the 
ratio of s to the pulling velocity v. In order to explain this uni- 
versal shape, we construct a simple geometric model for the 
interaction dynamics based on the observed stress axes. This 
simple model provides a quantitative prediction which repro- 
duces both the shape of the lens and its aspect ratio. 

Experiment: We perform experiments on sheets of vari- 
ous materials in a motorized frame which pulls one edge of 
the sheet away from the opposite edge. To start, we cut a pair 
of offset notches separated by a distance s\ these two cuts in- 
stigate the cracks which propagate inward and interact under 
uniaxial strain until they cut out a lenticular fragment, visi- 
ble in Fig. [If. The sheets are approximately 10 x 20 cm^ 
and pulling velocities are 0.7 to 14 mm/s. We collect data on 
the crack propagation and interaction using three techniques: 
shape analysis of the final fragment, crack-tip tracking to cap- 
ture the dynamics, and photoelastic imaging of the stress axes. 
For the final shape analysis, we remove the fragment from the 
apparatus and perform measurements on an unstressed sam- 
ple. We track the location of the tip on a sequence of images 
taken every 4 s (see Fig.[T])-f); during this tracking, the sam- 
ple remains in a stressed state. For the gelatin samples, we 
can use a polariscope to visualize the internal stresses during 
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FIG. 2: [Color Online] (a) Measurement of aspect ratio F of final 
fragment shapes for constant-v fracture of gelatin, cork, polystyrene 
foam, and aluminum foil, and uncontrolled fracture of nitrile. Dotted 
line is L = 2FF. (b) Average of edges i(w) for 8 gelatin fragments, 
with both lengths scaled by s. Dashed line is EcL[T]with A — \ 



and a = 2' which corresponds to the model (Eq.Mj). Inset: Sample 
image of fragment showing length L, width W '^ s, and extracted 
edge li(w). (c) Power law exponent a as a function of FF, obtained 
from fitting Eq.fl] dashed line is mean (a) = 0.49. (d) Power law 
prefactor A as a function of FF, obtained from fitting Eq.[T]with fixed 
a — \\ dashed line is mean (A) = 0.97. In (c,d) the paired points 
connected by vertical bars correspond to the two sides of a single 
fragment. All fits to Eq.fl] are performed over the portion oi li(yS) 
with wj s > 0.1 



fracture f26l, as shown in Fig.[T^h. 

The experiments investigate the effect of the initial crack 
offset 5, thickness h, material properties, and pulling veloc- 
ity V. Our primary material is 7%w gelatin {E ^ 10^ Pa, 
V = 0.4 to 0.5, h = 8 mm) 1271 . with comparison runs in 
nitrile (E - 10^ Pa, u - 0.5, h = 0.1 mm) |28|, cork 
(E ~ 10^ Pa, i^ ~ 0, h = 3 or 6 mm), polystyrene foam 
(E - 10^ Pa, i/ - 0, /i = 6 mm) |29|, and aluminum foil 
(E - 10^^ Pa, u - 0.3, h = 0.02 mm). These materials are 
selected for their diversity of material properties, and cover 
the full range of Poisson ratio values from to 0.5. For the 
gelatin, which is weak enough to break under its own weight, 
we support the sheet from below using a clear sheet lubricated 
with vegetable oil or water to minimize the frictional inter- 
action with the substrate. The ends of the gelatin sheets are 
affixed to the apparatus by means of mesh brackets cast into 
the sheet during preparation. All other materials are clamped 
at their ends and hang under their own weight; as a result, 
the sheet is not constrained to remain in a single plane during 
fracture. 



Shape: For each final fragment, we measure the length 
L and width FF, with FF ~ s corresponding to the initially- 
parallel crack edges (see Fig. [TJT). As shown in Fig. [2^, we 
observe an aspect ratio F = L/W ^ 2, independent of s, 
material, and v. Remarkably, F is constant even for runs in 
which the sheet curved out of plane, started from non-parallel 
initial cracks, or was subjected to shear in addition to the ten- 
sile strain. For the last two cases, the cracks reoriented to 
two straight paths perpendicular to the direction of loading, 
thereby still resulting a fragment with F ^ 2. 

To examine the universality of the shape which resulted in 
F = 2, we examine 8 gelatin fragments scanned at 600 dpi 
resolution. For each, we extract the coordinates [^, w] along 
the edge connecting the midpoint of the fragment to the tip, as 
shown in Fig. [2J3. Note that the two initial cracks (separated 
by s) remain approximately parallel until they reach this mid- 
point, after which they curve towards each other along the path 
i{w). These dynamics will be further quantified below. We 
examine the shape of the fragment by scaling both coordinates 
of £{w) by s and averaging over the ensemble of fragments. 
Each half-fragment spans < i^ < 1 and < f < 1; the 
full length is 2£, corresponding to F = 2. The average shape 
{£{w)) has an approximately power-law shape over much of 
its length, with the exponent varying from | initially to ^ for 
w/s>OA. 

Motivated by this average shape, we fit each edge (two from 
each fragment) to the functional form 
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for the region w/s > 0.1 and find mean values of (A) = 
0.97 ± 0.03 and (a) = 0.49 ± 0.01, with no systematic de- 
pendence on 5, material properties, or v (see Fig. |2]:d). In- 
dividual fragments vary from this mean behavior, particularly 
when the initial cracks were imperfectly parallel or the gelatin 
sheet was poorly lubricated. Thus, we conclude that the initial 
crack offset s sets not only L and FF of the final fragment, but 
is the only input to the universal shape function given in Eq.[T] 
Dynamics: To characterize the dynamics which result in 
this characteristic shape, we track the two crack tip positions 
for runs started from different s and pulled at three different 
V. We record the coordinates [xi{t)^yi{t)], where z = 1 for 
the stationary side and i = 2 for the pulled side. We set the 
origin of our moving coordinate system so that [x^y] = at 
the midpoint between the two tips, and take t = at the in- 
terpolated time when the two crack paths cross ^ = 0. The 
offset between the crack tips is measured at t = to be 
Axc = :^^i(0) — ^2(0) « s, which can be influenced by slight 
repulsion near t = (161 [20l. In addition, the characteristic 
timescale r = Axc/v is used to make comparisons between 
runs at different v, and we scale all distances by Axc, The 
resulting scaled variables are 
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FIG. 3: (a) Xi{i) and (b) yi(i) for 12 experiments on gelatin at '^ = 
0.7, 1.4 and 2.9 mm/s. (c) Pre-interaction stage: yi{t) for a single- 
crack (v = 0.7 mm/s) and a double-crack (t' = 1.4 mm/s). ^o < is 
the starting coordinate of each crack at time to < 0. (d) Interaction 
stage: {y{i)). 



In tracking the dynamics of the tip, all measurements are made 
on a stressed sample. 

Fig. [3^b examines both individual (thin lines) and aver- 
age (thick lines) trajectories, demonstrating that the scaling in 
Eq. [2] collapses the curves into a universal shape; no system- 
atic dependence on Axc or v is present in the scatter. Note 
that x{i) is on average constant for t < 0, indicating that little 
crack curvature is taking place at this stage. This results in 
each final fragment containing two parallel edges at opposite 
ends of the fragment. In addition, we observe that each crack 
initially behaves as if it were a single, isolated crack pulled at 
half the velocity, as shown in Fig. [5]:. For either a single crack 
or one of a crack pair, y oc {t -\-to)^^^. Therefore, we refer to 
t < as the pre-interaction stage and t > as the interaction 
stage. 

To examine how the observed similarity in £{w) arises 
from the dynamics, we consider {xi{i)) and {yi{i)), where 
(•) is the average over the 12 examples. During the interac- 
tion stage, (y) is observed to grow only logarithmically (see 
Fig. [3]l). This slowing down can be understood in terms of 
the rotation of the fragment, rather than fracture, becoming 
the dominant means of accommodating strain the thinner the 
remaining connection becomes (see Fig.[T]F.) 

These dynamics highlight several important features of the 
en passant crack geometry which may be quite general. First, 
the interaction between the cracks starts only when the tips 
pass each other. From that point, all length scales are set by 
the crack offset s and all time scales are set by the ratio of 
this offset to v. The final shape has an aspect ratio F = 2 
for materials with diverse elastic properties, including elastic 
modulus, Poisson ratio and heterogeneity. The curved portion 
of the fragment takes a power law shape with an exponent 
close to ^ . 

This universal shape suggests that a geometric model can 
provide insight into the dynamics, without reference to the 
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FIG. 4: [Color Online] (a) Model solution (thick lines) during inter- 
action stage (i > 0), with definitions of variables. Model in (b) final 
state and (c) translated from final state to show crack edges super- 
imposed on experimental image. Comparison to A = 1,3 shown as 
thin lines in (c). 



full elastic problem. Initially, the crack paths remain straight 
and parallel because there is a large bulk region between the 
cracks, allowing the principal axis of the stress to be aligned 
with the boundary loading. This is pure mode I fracture, and 
proceeds identically to a single crack, without torque (see 
Fig. [3]:). However, once the crack tips pass each other, the 
principal axis of the stress within the central region rotates 
and connects the tips. This provides a net torque on the cen- 
tral part of the sample and leads to locally mixed mode I and II 
loading. At each strain increment, the line connecting the two 
tips rotates and the subsequent fracture occurs relative to this 
new line. These features can be observed in the polariscope 
images of the internal stresses show in Fig.[T^h. 

Model: We construct a simple geometric model of the 
two-crack interaction for t > starting from three assump- 
tions. First, that the trace of the stress tensor is is largest along 
the line connecting the two cracks. Second, that each crack 
propagates orthogonal to this direction. This corresponds to 
the principle of local symmetry, whereby mode I is selected 
over mode II (301 [SH. These two assumptions alone are suffi- 
cient to generate curvature in the crack tip trajectory: once the 
cracks have moved past each other, the line connecting the tips 
rotates, and the crack direction rotates with it. For simplicity, 
we also assume that the two tips remain laterally separated by 
their initial offset s. As can be seen from Fig.|3^, x remains 
in the vicinity of ^ for t > 0, but with additional dynamics 
which we do not capture with the simple model. 

To formulate this model algebraically, we establish a co- 
ordinate system (see Fig. |4^) around the center of symmetry. 
Within this frame, the two edges are pulled away at equal and 
opposite velocities. During each strain interval, each outer 
crack edge (Rq) moves a distance dx away from the center. 
By continuity, each inner edge (Ri) moves in the opposite di- 
rection by the same amount. Each crack propagates forward 
to relieve the resulting stress, with the direction set by the or- 
thogonality condition and a distance chosen to keep each tip 
at fixed X— coordinate ±|. At each strain interval, the coor- 



dinates of the left-side tip are r = [— f , ^] and the cracking 
direction is normal to this vector, in the direction ft = [1, ^]. 
Thus, for each outward displacement by dx, the outer edge 
Ro = [x^ 0] will move leftward by —dx so that x ^ x — dx, 
the inner edge Ri = [— (5 + x), 0] will move rightward by dx, 
and the y coordinate of the tip will advance by y ^ y -\- ^ 
to correspondingly relieve the strain. This provides the differ- 
ential equation 



d^ ^ ^ 
dx 2y' 



(3) 



Integration yields (|)^ = - 
integration is set by y( — |) = 
imental results, we define y 
leads to the simple shape 

y 



(I + f ), where the constant of 
= 0. To compare with the exper- 

= ^ and X = -(^ + ^), which 



: X2. 



(4) 



The crack propagation stops when the inner crack edge meets 
the tip of the other crack, which sits at | . This occurs when 
— {s-\-Xf) = f , or equivalently Xf = yf = 1. Impor- 
tantly, this solution matches the F = 2 shown in Fig. [2^ and 
(A) = 0.97 and (a) = 0.49 shown in Fig. ^d. For com- 
parison, we have superimposed the model prediction on an 
experimentally-obtained image in Fig. [4]: and find excellent 
agreement except at the early stages of curvature. This early 
regime corresponds to the portion of the shape with a = | in 
Fig.[2p, and to the time immediately after i = Fig.jSp. 

While we have constructed the model in terms of boundary- 
driven strain increments dx, Eq. [3] and the same final shape 
would result for any brittle fracture mechanism, including 
cyclic thermal expansion and contraction, as long as the crack 
growth still occurred in the direction normal to the line sep- 
arating two tips which remained at fixed separation. For a 
boundary loading which provides a different stress intensity 
factor ratio than observed here, the crack may grow at a larger 
or smaller (non-orthogonal) angle with respect to r |[T6ll . In 
such a case, the factor 2 in Eq. |3] can be generalized to a 
constant A with —^ = j-. When integrated, the solutions 

for A 7^ 2 retain a = ^ but have A = A{X) and therefore 
r 7^ 2, as shown in Fig. [4]:. This last case is noteworthy, since 
in geologic observations, F > 2 is commonly observed |4|. 
In experiments where we provide additional lateral boundary 
stresses, we are able to achieve both larger and smaller as- 
pect ratios. This raises the possibility of using lenticular crack 
shapes as a diagnostic for the stress conditions under which 
cracks were formed in nature |32|. In particular, F may serve 
as a means to infer the boundary loading in situations where 
history and dynamics are inaccessible, for instance in obser- 
vations of Europa's ice sheets O. 

In conclusion, we have identified the length and time scales 
which control a generic two-crack interaction problem. A 
simple geometric model allows us to identify the key mech- 
anism leading to the universal shape: the largest tensile stress 



follows a straight line connecting the tips and the cracks prop- 
agate perpendicular to this line. 
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